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ABSTRACT 

We show that recent observations of He I and N II lines of 77 Carinae may provide support for an 
orbital orientation where the secondary star is closest to us at periastron passages. This conclusion 
is valid both for the commonly assumed masses of the two stars, and for the higher stellar masses 
model where the very massive evolved primary star mass is Mi = 170M 0 and its hot secondary star 
mass is M 2 = 80 Mq. The later model better explains the change in the orbital period assuming that 
the ninetieth century Great Eruption was powered by accretion onto the secondary star. Adopting 
the commonly used high eccentricity e — 0.9 and inclination i = 41°, we obtain a good fit to newly 
released Doppler shift observations of He I emission and absorption lines assuming they are emitted 
and absorbed in the acceleration zone of the secondary stellar wind. Our model in which the secondary 
star is in the foreground at periastron is opposite to the view presented recently in the literature. 
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1 . INTRODUCTION 

77 Car is a binary system (Damineli 1996; 
Damineli et al. 1997) composed of a very massive 
primary star (Davidson & Humphreys 1997) and 
a hotter and less luminous evolved main sequence 
secondary star. Despite two decades of detailed obser¬ 
vations (e.g., Smith ct al. 2000; Duncan & White 2003; 
Whitelock et al. 2004; Corcoran 2005; Davidson et al. 
2005, 2015; Smith 2006; Hamaguchi ct al. 2007; 

Damineli et al. 2008b; Corcoran et al. 2010; 
Martin et al. 2010; Mehner et al. 2010, 2012, 2015; 
Abraham et al. 2014; Hamaguchi et al. 2014a,b ) 
and modelling (e.g., Pittard et al. 1998; Soker 2001; 
Pittard & Corcoran 2002; Akashi ct al. 2006, 2013; 
Kashi & Soker 2009a; Okazaki ct al. 2008; Smith 
2010; Groh et al. 2012; Madura et al. 2013, 2015; 
Clementel et al. 2015b) there arc disagreements over 
two important properties of the binary system. These 
are the masses of the two stars and the orientation of 
the eccentric orbit. 
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The low-masses model assumes that rj Car is at its Ed¬ 
dington luminosity limit and a mass of > 120 Mq is de¬ 
rived (Hillier et al. 2001). Some studies take it to be the 
mass of the primary star, while other assumes that the 
combined masses of the two stars amounts to that mass 
with Mi > 90Mq and M 2 > 30 Mq (e.g., Okazaki et al. 
2008; Clementel et al. 2015b). 

The high-masses model was developed by us 
(Kashi & Soker 2010a) under the assumption that most 
of the extra energy released during the 1837-1856 Great 
Eruption (GE) of 77 Car originated from high accretion 
rate onto the secondary star. There are models that 
explain the GE with only one star (e.g., Matt & Balick 
2004), and others that use three stars (Livio & Pringle 
1998; Portegies Zwart & van den Heuvel 2016). 

The GE was agiant eruption and a SN Impostor, 
which is the group of eruptive massive stars within 
the more general groups of intermediate luminosity 
optical transients (ILOTs; Kashi & Soker 2015). In 
Kashi & Soker (2015) we presented the high-accretion- 
powered ILOTs (HAPI) model, according to which all 
ILOTs are powered by high accretion rate onto a main 
sequence (MS), or a star slightly evolved off the MS 
(Kashi & Soker 2010b). According to the HAPI model, 
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then, the luminosity peaks of the Great Eruption re¬ 
sulted from accretion onto the companion close to peri- 
astron passages. Those peaks are « 5.1 5.2 years apart. 
If one assumes that the peaks are related to peristron 
passage, and occur about the same time after (or be¬ 
fore) periastron, the separation between the peaks can 
be inferred as the orbital period. The second peak came 
after the eruption has started, and mass was lost. There¬ 
fore the orbital period is smaller than the time interval 
between the peaks. We adopted a value of 5.1 years. 
The present orbital period is 5.54 years, suggesting that 
during the Great Eruption the orbital period has in¬ 
creased. The orbital period changed as the secondary 
accreted mass, and as mass was lost from the binary 
system, both as a wind from the primary star and as 
jets from the secondary star. 

In that study (Kashi & Soker 2010a) we found that 
for the HAPI mechanism to work within the orbital 
changes constraints, the two stars should have a signif¬ 
icantly larger masses than 12OM 0 . The masses should 
be in the range of Mi ~ 15O-2OOM 0 and M 2 — 60- 
9OM 0 . This high-masses model is further supported 
by evolution of massive stars on the HR diagram. The 
calculations of Figer ct al. (1998) show that a zero-age 
main sequence (ZAMS) star with an initial mass of 
A?zams — 23OM 0 is required to explain the present 
luminosity of the primary of 77 Car. Indeed, modern 
stellar evolution tracks of 12OM 0 stars even when ro¬ 
tation is considered (Ekstrom et al. 2012; Georgy et al. 
2012 ) do not reach that luminosity, and more massive 
models of Mzams — 250M 0 are required, as found by 
Chen et al. (2015). The primary lost large amounts of 
mass to the bipolar nebula around 77 Car, the Homuncu¬ 
lus (Smith & Ferland 2007; Gomez et al. 2010). When 
the large mass loss along the evolution is considered, a 
present mass of ~ 150 — 2OOM 0 is compatible with the 
ZAMS mass inferred above. 

In the debate on the orientation of the binary 
system, the two sides are holding a literally ~ 180° 
opposite views. One side holds that during periastron 
passages the primary star is closer to us (w ~ 240°- 
270°; e.g., Iping et al. 2005; Nielsen et al. 2007; 

Damineli et al. 2008b; Henley et al. 2008; Parkin et al. 
2009; Groh et al. 2010; Gull et al. 2011; Madura et al. 
2012; Clementel et al. 2015a; Richardson et al. 2015; 
Teodoro et al. 2016, while the other side holds the 
view that during periastron passages the secondary 
star is at its closest location to us (w ~ 90°; e.g., 
Abraham et al. 2005; Falceta-Gongalves et al. 2005; 
Abraham & Falceta-Gongalves 2007; Kashi & Soker 
2008, 2009b, 2011; Tscbrenko et al. 2013). 

The disagreement on the orientation stems mainly 
from the controversy on the source of some emission 
and absorption lines, in particular the He I lines. In 


a recent paper Richardson et al. (2015) present a study 
of the variations of some spectral lines with the phase 
of the binary orbit, based on observations with the 
CTIO 1.5m telescope. They attributed the He I lines 
to the primary star. In Kashi & Soker (2007), on the 
other hand, we assumed that the He I lines observed by 
Nielsen et al. (2007) are formed in the acceleration zone 
of the wind blown by the secondary star. We calculated 
the Doppler shift variations of the lines as a function of 
orbital phase with the low-masses model of 77 Car, and 
found that a good fit is obtained if the He I lines are 
formed in the region where the secondary wind speed is 
Vzone = 430 km s- 1 . 

In the present study we combine the high-masses 
model of 77 Car with the assumption that the He I lines 
originate in the acceleration zone of the secondary stel¬ 
lar wind, and try to fit the new Doppler shifts presented 

by Richardson et al. (2015). 

2 . THE CONTROVERSY OF THE HELIUM I LINES 

The origin of the visible He I P Cyg lines (A7065A, 
A5876A, A5015A, A4992A, and A447lA) in the binary 
system 77 Car is in dispute, with different researchers at¬ 
tributing it to different regions, e.g., the primary star 
(e.g., Falceta-Gongalves et al. 2007; Humphreys et al. 
2008; Richardson et al. 2015). We attribute the He I 
P Cyg lines to the acceleration zone of the secondary’s 
wind (Kashi & Soker 2007, 2008). A similar dispute 
exists for the N II AA5668-5712A line (Mehner et al. 
2011a). While Mehner et al. (2011a) argued that it can¬ 
not come from the secondary star, we noted that the 
N II lines closely follow the behavior of the He I lines, 
and attributed it to the secondary wind ( ishi k Soker 
2011 ). 

The effective temperature of the secondary was es¬ 
timated be T e ff 2 ~ 34 000-38 000 K by Verner et al. 
(2005), and more recently T e g i2 ~ 40 000 by 
Mehner ct al. (2010). The main arguments for the 
He I P Cyg lines origin in the acceleration zone of 
the wind blown by the secondary star are as follows 
(Kashi & Soker 2011). 

1. The Doppler shift of the P Cyg absorption com¬ 
ponents follows very well the secondary’s orbit, as 
we show in section 3 for the high-masses model, 
and as was shown before for the low-masses model 

(Kashi & Soker 2007, 2008). 

2. The lines are known to originate in stars with tem¬ 
peratures well above 30 000 K, mostly in hydro¬ 
gen deficient stars (e.g., Leuenhagen et al. 1996; 
Crowther & Bohannan 1997; Grunhut et al. 2013; 
Wessolowski et al. 1988; Leuenhagen & Hanrann 
1998; further discussion is in section 5). 
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3. The secondary can account for the amount of ab¬ 
sorption in the He I lines, as we show in section 

4. 

4. The Doppler shift of the emission follows that of 
the absorption, as we show in section 3. 

The secondary mass loss rate is much lower than the 
primary’s. It is important to mention that most of the 
lines observed from the 77 Car system do originate in the 
primary and its wind (the best example is probably the 
hydrogen lines; Weis et al. 2005, Davidson et al. 2005). 
Take H I A4103A line for example. This line, and many 
others, originate in the primary. But it shows a P Cyg 
profile that shifts much less in periastron compared to 
the He I P Cyg lines we discuss here. 

It might seem a “strange” coincidence that the place 
in the secondary wind where the He I are formed, ac¬ 
cording to our model, has about the same velocity as the 
primary’s wind. However, even a stranger coincidence 
exists for a model where the lines are formed in the pri¬ 
mary’s wind: the area that absorbs the lines would have 
to change its velocity in the same way the secondary 
moves around the center of mass of the binary system. 

Using the assumption that the He I lines originate in 
the secondary wind, in the past we fitted their Doppler 
shift variations with orbital phase for the low-masses 
model (the conventional model) of 77 Car, with Mi = 
120 Mq and M 2 = 30M 0 (Kashi & Soker 2008). The 
Doppler shift of the P Cyg absorption component of the 
He I lines was found to be in agreement with the binary 
orientation with a longitude angle oj = 90°, i.e., sec¬ 
ondary closest to us at periastron. We here upgrade the 
model to include the new Doppler shifts presented re¬ 
cently by Richardson et al. (2015), and the high-masses 
model that better fits the luminosity and the behavior of 
77 Car during the Great Eruption according to the HAPI 
model. 

In fitting the He I lines Doppler shifts with orbital 
phase we scan the following parameter space. 

1. The eccentricity is e ~ 0.85-0.93. 

2. Inclination angle (the angle between a line perpen¬ 
dicular to the orbital plane and the line of sight) 
is also in the consensus to be * cs 41°. 

3. The orbital period is P = 2023 days 

(Damineli et al. 2008a). 

4. For the masses we use our results (Kashi fc Soker 
2010a) that the present masses are Mi ~ 150- 
200M e and M 2 — 60-90 Mq. We take here for the 
present masses of 77 Car components Mi = 17OM 0 
and M 2 = 80 Mq. 


5. We assume that the observer is behind the sec¬ 
ondary at periastron, namely w ~ 90°. 

We note that there is inconsistency in some papers 
regarding the exact epoch of periastron, or phase 0 . 
Richardson et al. (2015) used JD 2454842.5 for the 2009 
periastron passage. Mehner et al. (2011a) used JD 
2454860 as a reference time for the event (see also dis¬ 
cussion at the appendix of Mehner et al. 2011b). The 
uncertainty in determining the time of periastron may 
cause the change in radial velocity to appear after pe¬ 
riastron rather than before, or vice versa. Evidently 
the observations collected so far from 77 Car are insuffi¬ 
cient for determining the exact time of periastron. We 
will adopt an intermediate value between the references 
above, i per = JD 2454850. This value coincides with 
having the sharp variation in radial velocity of the lines 
studied here at periastron. However, one should bear 
in mind that the uncertainty of the periastron epoch is 
of ±10 days. We corrected the orbital phases inferred 
from the observation dates in this work to have phase 0 
at t per- 


3. ORBITAL PARAMETERS 


Based only on geometric considerations, i.e., neglect¬ 
ing variations of the wind speed near periastron passage 
and stochastic wind speed variations, we here use our 
model to fit the observations of Richardson et al. (2015). 
The orbital velocity of the secondary relative to the pri¬ 
mary, v OI b , is converted to the velocity relative to the 
center of mass 
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where f(t) = r(t)/a. 

The entire stellar mass dependency is embedded in 
the factor For the low-masses model (Mi = 120 Mq 
and M 2 = 30 Mq) its value is /m = 4.25M)/ 3 , while for 
the high-masses model (Mi = 17OM 0 and M 2 = 80 Mq) 
it is f m = 4.28 Mq 3 . Therefore the amplitude of the fit 
with the new parameters is less than 1 % larger than the 
amplitude using the old parameters. 

On top of v m there are factors related to the obser¬ 
vation angle w, the inclination i. and the eccentricity e. 
From the result we subtract the constant velocity 77 ZO ne 
of the zone in the secondary wind where the lines are ab¬ 
sorbed, or the average velocity of the emitting gas when 
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the peak emission is fitted. Figure shows in a blue-solid 
line our fit to the observed radial velocity absorption 
component of He I lines from Richardson et al. (2015) 
and of Nielsen et al. (2007). It is clear from Figure 
that our fit to the new lines at A4922A and A5015A is not 
as good as our fit to the lines studied by Nielsen et al. 
(2007). We attribute this to contamination of the two 
A5015A and A4922A lines by Fe II lines. This contami¬ 
nation was mentioned by Richardson et al (2015). 

More observations of He I lines, specifically the 
He I A4714A line, were taken for the 2009 event by 
Mehner et al. (2011b), and for the 2014.6 event by 
Mehner et al. (2015). The later paper summarizes ob¬ 
servations of that line from the previous 3 events (2003.5, 
2009, and 2014.6). Figure 2 shows how our model fits the 
observations. We used the same model we used for the 
other He I lines in Figure 1, but u zon e,abs = 370 km s -1 . 
We take this slightly different value than the value of 
fzone.abs = 430 km s _1 that was used in Figure 1 to 
match the average of the absorption component’s ra¬ 
dial velocity. The different values mean either that 
the He I A4714A line is absorbed in the wind slightly 
closer to the wind origin on the secondary star, or that 
there are large variations and uncertainties in the de¬ 
rived Doppler shifts near periastron passages. 

Mehner et al. (2015) found that the 2014.6 event was 
different than pervious events, and according to their 
interpretation in the framework of the accretion model, 
it showed signs of less accretion onto the secondary close 
to periastron, indicating weaker primary wind. The He I 
line flux increased significantly in 2009-2014 compared 
to 1998-2003. The absorption of the He I A4714A line 
disappeared 8 days before periastron and re-appeared 8 
days after (assuming our t per ). The radial velocities are 
lower in the 2014.6 event compared to previous events. 
In fact, for fitting only the data from 2014.6 it would be 
better to use fzone.abs = 330 km s , keeping the rest 
of the parameters unchanged. This behavior, together 
with the large fluctuations in the Doppler shifts near pe¬ 
riastron as seen in the different figures, suggest that the 
velocities of the regions where lines are formed in the dif¬ 
ferent regions of the secondary wind vary from cycle to 
cycle and in short time scales near periastron passages. 
We should therefore aim at fitting the general behavior 
of the Doppler shift variations with orbital phase; a per¬ 
fect fit to the Doppler shifts of this interacting binary 
system cannot be achieved using the geometrical effects 
alone. 

We used the same principles we used for fitting the ra¬ 
dial velocity variations of the He I absorption lines to fit 
the emission peak bisector velocity of the He I A6678A 
line across the 2009 event observed by Richardson et al. 
(2015). We find that a value of » ZO ne,emi = 60 km s -1 
for the average velocity of the gas emitting the line 



orbital phase 



Figure 1. Our fit to the observed radial velocity absorp¬ 
tion component of He I lines. The black diamond points are 
the He I A7065A, A5876A, A5015A, A4992A, and A447lA 
lines from n al (2( ) with the original error bars. 

The CTIO radial velocity data from >n (2( ) 

are shown in green circles (He I A4922A) and red squares 
(He I A5015A), and include no error bars. We added bars to 
two points of each set that indicate our estimate of the fluc¬ 
tuations in the velocity of the wind absorption component, 
±40 km s -1 . Observation phase was modified to periastron 
time defined in the text. Our fiducial model assumes that 
the spectral lines originate in the secondary stellar wind, that 
to = 90° (i.e., secondary star closest to us at periastron), an 
eccentricity of e = 0.9, and stellar masses of Mi = 170Mq 
and M 2 = 80 Mg. The upper panel shows the entire orbit, 
and the bottom panel is zoomed to times near periastron. 
The lower panel includes also a model with e = 0.93 instead 
of e = 0.9 (grey line). For comparison we present a fit for a 
case where the He I lines originate from the primary stellar 
wind, and the orientation is to = 270° (dashed-dotted pink 
line). The velocity of the region in the wind responsible for 
the lines is taken to be u ZO ne,abs = 430 km s -1 . 

gives the best results. Figure 3 shows the fit we ob¬ 
tained, together with the observations of the emission 
peak bisector velocity of the He I A6678A. We present 
two other cases with cu = 90° (secondary star closer to 
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Figure 2. Like Figure , but taking u zone , a bs = 370 km s -1 
to fit the observed radial velocity absorption component of 
the He I A4714A line, from observations of 

(2011b) and Mehner et al. (2015). 


the observer at periastron), one with an eccentricity of 
e = 0.93, and the other is the low-masses model from 
Kashi & Soker (2008). We also show two models where 
the line is assumed to originate in the primary stellar 
wind. One is our simple geometric model, but the line 
is emitted by the primary stellar wind and w = 270° 
(dashed-dotted line), and the second is the fit presented 
by Richardson et al. (2015). It is clear that the models 
where the line is emitted by the secondary stellar wind 
and the secondary star is in the foreground at periastron 
result in a much better fit to the observations. 

We next test our model against the variation 
in the Doppler shift of the absorption component 
of the N II A5666-5771A lines, taken from both 
Richardson et al. (2015) and Mehner et al. (2011a). We 
use here u ZO ne,abs = 370 km s -1 . The N II A5668A 
line observed by Mehner et al. (2011a) has clear emis¬ 
sion peaks close to periastron, but these were not clear 
in the profiles of Richardson et al. (2015). We fit also 
the emission component of the N II A5668A line, with 
Dzone.emi = 70 km s . The fits are presented in Figure 
4. 

Though the N II and He I lines intensities behave dif¬ 
ferently, the radial velocity of the N II A5668A emis¬ 
sion component follows that of the absorption compo¬ 
nent amazingly well. Also, the entire lines-formation 
regions show the same Doppler shift variation. Such a 
behavior cannot be explained in frame of a model where 
the N II and He I lines-formation regions change their 
location within the primary wind. 

The main driver of the variations in the Doppler shifts 
is the pure orbital motion. However, it is not the only 
one. A variation in location within the secondary wind 
in our model can take place, and it is even expected to 


Figure 3. Our fit to the emission peak bisector velocity 
of the He I A6678A line (blue line). We show also the low- 
masses model (yellow line), and a higher eccentricity model 
(grey line). Bars indicate our estimate of fluctuations in 
the average velocity of the emitting gas in the secondary 
wind. Also shown is the case where the line originates in the 
primary stellar wind and the primary star is closest to us at 
periastron (dashed-dotted pink line). The black dashed line 
is the fit to the data taken from ;ou t (2015), 

who assume that the line originates in the primary wind. 


occur. As the secondary approached periastron, its wind 
properties change, most likely due to accretion of gas 
from the primary stellar wind (Soker 2005; Akashi et al. 
2013). The accretion phase lasts for several weeks, 
but the influence on the secondary stellar wind prop¬ 
erties can last few months after the event, until the ac¬ 
creted mass is removed by the restoring secondary wind 
(Kashi & Soker 2009). It is therefore very reasonable 
that lines will be absorbed in different locations across 
periastron. This secondary effect is the main cause for 
the deviation from our Doppler-based model, as seen in 
the figures. 


4. ABSORPTION OF HELIUM BY THE 
SECONDARY 


According to Nielsen et al. (2007) the amount of ab¬ 
sorption of the He I A5015A P Cyg line reaches up to 
50%. This line, however, shows different orbital change 
(probably because it is blended) and is not one of the 
lines we claim to originate in the secondary wind (see 
section 3). The He I lines we attribute to the secondary, 
show much less absorption. The He I A7067A P Cyg 
line, for example, shows 10% absorption near periastron 
Nielsen et a' (2007). We can use this to better constrain 
the properties of the two stars. 

The maximum absorption, assuming it occurs in the 
secondary wind, is obtained when 


A 


max 


A2-L2 

fhLl + / 22 -L 2 


( 3 ) 
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Figure 4. Upper panel: Our fit to the observed ra¬ 
dial velocity absorption component of the N II A5666A and 
N II A577lA lines from the 2009 event (Richardson et al. 
20 ), and the N II A5668A line observed by M< ner al 

(201. ) in the 3 previous events. Lower panel: Our fit 
to the observed radial velocity emission component of the 
N II A5668A line from Mehner et al. (2011a). 

Where L\ and L 2 are the luminosities of the primary 
and the secondary, respectively, and fn and // 1 are their 
fractions within the line absorption waveband. 

According to Nielsen et al. (2007), the He I A7067A 
shows 10% absorption for a range between —700 
—500 km s -1 . Taking A max < 0.1 in equation (3) we 
get the requirement that 


Assuming Blackbody radiation, and taking conventional 
parameters for the temperatures of the primary and the 
secondary Ti = 20 000K and T 2 = 40 000K, respectively, 
we obtain fn/f 12 = 5.9. This means that according 
to equation (4), the luminosity ratio needs to satisfy 
L 2 /L 1 > 0.65. According to the conventional parame¬ 
ters, the ratio L 2 /L\ = 0.2, so this leads to a contra¬ 


diction. It is not possible to find for the conventional 
model, a point in the stellar evolution path, even with a 
different effective temperature, that would be even close 
to satisfy that requirement. We therefore conclude that 
the two stars must have a smaller luminosity ratio, and 
consequently smaller mass ratio, as is expected in the 
massive stars model. They should also have a smaller 
temperature ratio to make the requirement in equation 
(4) easier to meet. 

For the massive stars model we propose, we can track 
the evolution of M\ = 170Mq and M 2 = 80 Mq stars 
(e.g., Ekstrom et al. 2012, Kohler et al. 2015 and refer¬ 
ences therein), and try to find a reasonable set of param¬ 
eters that satisfies the above equation. We find that if 
we take T x = 25 000K (Hillicr et al. 2001), T 2 = 37000K 
(Verner et al. 2005), which gives fn/fi 2 = 2.8, then the 
requirement from equation (4) becomes L 2 /Li > 0.31. 
And indeed, stellar evolution tracks give L\ ~ 3 x 10 6 Tq 
and L 2 — 1.2 x 10 6 Lq, satisfying the condition. 

5. SUMMARY AND DISCUSSION 

We used spectroscopic observations of 77 Car close 
to the 2009 periastron passage (Richardson et al. 2015; 
Mehner et al. 2011a,2011b) to show they support ear¬ 
lier suggestions that companion is in the foreground at 
periastron ( vaslii k Soker 2011 and references therein). 

We assumed that the He I and N II spectral lines 
originate in the acceleration zone of the secondary star 
of rj Car (Kashi & Soker 2007, 2011). We then took 
the secondary to be closest to us at periastron pas¬ 
sages, i.e., uj = 90°. We further used the high-masses 
model of 77 Car with component masses of Mi = 170M Q 
and M 2 = 80 Mq. These masses better fit evolutionary 
tracks of massive stars that cross the locations of the 
two stars on the HR diagram, than the commonly used 
masses of Mi = 120Mq and M 2 = 30Mq. A massive 
primary star with initial mass larger than Mi = 200Mq, 
and a secondary star with an initial mass larger than 
M 2 = 50Mq are supported by stellar evolution calcula¬ 
tions for very massive stars (e.g., Yungelson et al. 2008; 
Brott et al. 2011; Yusof et al. 2013). 

The high masses, more generally in the range of Mi cs 
15O-2OOM 0 and M 2 — 60-90 Mq, can account also for 
the powering of the nineteenth century Great Eruption 
by mass accretion onto the secondary star (the HAPI 
model; Kashi & Soker 2010a) . For the eccentricity and 
inclination angle we used the commonly accepted val¬ 
ues of e ~ 0.9 and i = 41°, respectively. As evident 
from Figs. 1-4 we could fit the general variation of the 
Doppler shifts with orbital phase. Therefore, the sug¬ 
gestion that the secondary star is in the foreground at 
periastron (e.g., Kashi & Soker 2008; Tsebrenko et al. 
2013) is definitely tenable. The suggestion that the pri¬ 
mary star is in the foreground at periastron seems to 
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encounter problems, e.g., as evident from the dashed 
black line in Figure 2 that is the model proposed by 

Richardson et al. (2015). 

In some cases the opposite model, of a line originat¬ 
ing from the primary stellar wind and the primary star 
is in the foreground at periastron passages, might fit 
part of the Doppler shifts, e.g., fitting the data from 
Richardson et al. (2015) in Figure 1. However, this 
model has too low amplitudes both at periastron pas¬ 
sages, as evident from all figures, and away from perias¬ 
tron passages, e.g., right side of Figure 2. 

In the present study we considered only the role of the 
orbital motion on the variation of the Doppler shift with 
orbital phase. It appears clear from the fluctuations in 
the Doppler shift values and from the non-perfect fit¬ 
ting that the velocity of the zone responsible for the 
formation of each line is changing. Both stochastic vari¬ 
ations and variation with orbital phase, noticeably near 
periastron passages, exist. These variations are another 
manifestation of the unrelaxed nature of this binary sys¬ 
tem. 

The three unknowns about the binary system dis¬ 
cussed here are: (i) masses of two stars; (ii) orbital ori¬ 
entation; (iii) exact periastron time. The masses of the 
two stars were obtained by Kashi & Soker (2010a) as 
explained in section 1. Here we showed that the mass 
estimate obtained by Kashi & Soker (2010a) also allows 
to fit radial velocities of the lines. The orbital orienta¬ 
tion is the main parameter discussed here in detail. The 
time of periastron is uncertain. We used a fixed value 
here and did not fine-tuned its value. It may however be 
possible to use the He I and other lines to get a better 
constrain for its time. What we find here is that our fits 
are consistent with periastron time of JD 2454850 but a 
few days difference is also possible. 

In addition to the Doppler shifts presented here, there 
are other arguments that support the suggestion that 
the secondary star is in the foreground at periastron 
passages. The four supporting arguments listed by 
Kashi X Soker (2008) include the evolution of the radio 
emission and the behavior of the He I A10830A line. In 
a previous study (Kashi & Soker 2009b) we further ar¬ 
gued that the column density toward the X-ray emitting 
gas, that is the post-shock secondary wind, is more com¬ 
patible with a binary orientation where for most of the 
time the secondary star is in the background, being on 
the foreground only near periastron passages. Another 
supporting argument was brought by Tscbrenko et al. 
(2013). They demonstrated that the asymmetric mor¬ 
phology of the blue and red-shifted components of the 
outflow at hundreds of astronomical units from 77 Car, 
can be accounted for from the collision of the free pri¬ 
mary stellar wind with the slowly expanding dense equa¬ 
torial gas closer to us. Namely, for most of the orbital 


period the primary is in the foreground, and at perias¬ 
tron passages the secondary star is in the foreground. 

Humphreys et al. (2008) discovered that the He I lines 
were absent from rj Car’s spectra prior to the mid 1940s. 
They discussed the various difficulties it may pose to 
models attributing the required He ionizing photons to 
the the secondary, ionizing the primary wind. They 
also conclude that even a much denser primary wind 
could not have obscured the He ionizing photons com¬ 
ing from the secondary. Clementel et al. (2015b) and 
Mehner et. al. (2015) claimed, however, that the He ion¬ 
izing photons only moderately penetrate the dense post- 
shocked primary wind. Therefore if the primary wind 
at that time was 2-4 times denser, it should have been 
enough to change the ionization structure of He in the 
primary’s wind, and prevent the formation of the lines. 
These arguments, however, are irrelevant for a model in 
which the He I lines originate in the secondary’s wind. 
In the frame of the accretion model, however, it is easier 
to provide an explanation. A denser primary wind can 
form a thick accretion belt around the secondary close to 
periastron (Kashi & Soker 2009) that would last for the 
entire orbit, providing a shield for its radiation. Even 
after the belt is gone, the mass that is accreted onto the 
secondary changes its photospheric structure and makes 
it cooler, diminishing the He ionizing photons. This 
may require that the primary wind before the 1940s was 
20-30 times denser. A new study by Kashi et al. (2016) 
suggests that the mass loss could have reached that mag¬ 
nitude or even higher. 

Mehner et al. (2011a) also observed the N II lines and 
their velocity shifts from reflected polar spectra at the lo¬ 
cation known as “FOS4”. They argue it may be a prob¬ 
lematic observation for the orbital motion explanation to 
the Doppler shift variations. A similar argument appear 
in Mehner et al. (2011b) regarding the He I A4714A line. 
The only direct comparison of FOS4 and direct view is in 
figure 8 of that paper. Though it may at first sight look 
like the radial velocity of the absorption in FOS4 follows 
that of the direct view, we notice the following: (1) The 
observations at —353 and —82 days (phase —0.18 and 
—0.045) show very high radial velocity in absorption, 
much above other observations of the same line. ( 2 ) At 
—82 days the value of FOS4 is 40 km s -1 higher. As the 
observations are sparse, there is now way to know if this 
is significant. It may indicate that the FOS4 Doppler 
shift of the absorption is smaller in amplitude and non 
systematic. Alternatively, it may be a fluctuation. It is 
therefore impossible to know if FOS4 persistently follows 
the direct view or not. A detailed comparison of densely 
sampled multiple lines is needed in order to check that. 
Even if it does, it may well be possible that not only po¬ 
lar light is reflected to FOS4 and there is some reflection 
from equatorial regions. 
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Mehner et al. (2010) used the distribution of gas and 
ionizing radiation around rj Car to constrain the prop¬ 
erties of the secondary. If the limits of Mehner et al. 
(2010) hold, then the secondary mass should be M 2 < 
60Mq. A mass of M 2 = 60M 0 is still within the HAPI 
model for the Great Eruption. We here showed that 
the Doppler shifts can be fitted with M 2 = 30M 0 and 
M 2 = 80Mq Any value for M 2 in this range can be 
fitted. But as stated, because of the luminosity of the 
primary star and the HAPI model for the GE, we prefer 
the high-masses model. 

P Cyg profiles in He I are found in hot hydrogen poor 
stars (Leuenhagen et al. 1996; Wessolowski et al. 1988; 
Leuenhagen & Hamann 1998). We know from observa¬ 
tions that the primary’s outer layers consists of about 
50% helium (Davidson et al. 1986; Dufour et al. 1997). 
A few Mq of material from the primary were accreted 
onto the secondary during the eruptions. The accreted 
gas makes the secondary’s envelope enriched with he¬ 
lium. It is very plausible that even though the secondary 
is hot, its helium lines are stronger than other stars in its 
evolutionary stage, dumpiness of the secondary wind 
(Kashi & Soker 2007) can make part of the gas some¬ 
what cooler, also enhancing the He I lines. 

The very massive primary star of Mi > 150M 0 and 
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the very high eccentricity of the binary orbit hints that 
the system was once a triple system, and that the pri¬ 
mary formed by the merger of two (or more) stars. The 
merger process released large amount of gravitational 
energy within weeks to months. Such an event can 
be classified as ILOT. Therefore, it may well be that 
the nineteenth century Great Eruption was not the first 
ILOT of this system. We note that the estimates for 
the masses in the HAPI model do not depend on the 
previous existence or nonexistence of a third star. The 
masses of only two stars are relevant for both the cal¬ 
culations done in this paper for spectral fitting, and the 
calculation done in Kashi & Soker (2010a) for modeling 
the light curve. The third star is suggested as a possi¬ 
ble easier route to obtain the large mass of the primary, 
together with a high eccentricity orbit. But it is well 
possible that both were obtained with only two stars 
along the entire evolution. 
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